Clinical Cancer Research

Purpose: Limited CNS bioavailability and pharmacodynamics are obstacles to effective systemic therapies for glioblastoma. One strategy to overcome these challenges is drug combinations enhancing CNS penetration and/or tumor chemosensitivity. LP-184, a synthetic acylfulvene class alkylator, induces DNA damage and inhibits glioblastoma cell viability in pre-clinical models. LP-184 is a prodrug converted to active metabolites by intracellular prostaglandin reductase 1 (PTGR1) that is over-expressed in >70% of glioblastoma. DNA damage induced by LP-184 is MGMT agnostic and reversed by transcription-dependent NER. Patients: LP-184 was evaluated in a Phase 1a study (NCT05933265) in 63 adult patients with advanced malignancies including 16 patients with recurrent glioblastoma. All patients with glioblastoma received prior standard-of-care therapy and most had received 1 or more additional therapies before enrollment. Results: Patients with glioblastoma experienced more frequent transaminitis, Grade 1-2 nausea and a trend towards more frequent and severe thrombocytopenia compared to the non-glioblastoma cohort. Otherwise, overall toxicity profiles were similar. Clinical pharmacokinetic analysis combined with published pre-clinical intra-tumoral bioavailability data (~20% penetration) predicted that LP-184 at the recommended dose for expansion (RDE) would achieve cytotoxic levels if combined with spironolactone, a BBB permeable ERCC3 degrader and TC-NER inhibitor that sensitizes glioblastoma cells to LP-184 3-6-fold. We show that three daily doses of spironolactone deplete orthotopic glioblastoma PDX ERCC3 protein by ~ 80% and increases tumor LP-184 cytotoxicity 2-fold. Conclusions: LP-184 is well tolerated at the RDE, and we establish a clinically translatable scheme for dosing spironolactone in combination with LP-184 for a future Phase 1b clinical trial.

1.The ABACUS study was a single arm, phase II trial evaluating neoadjuvant atezolizumab in operable urothelial carcinoma (UC). Initial bulk transcriptomic and immunohistochemistry analyses suggested links between immune activation, tissue remodeling, and resistance pathways such as transforming growth factor {beta} (TGF{beta}) that were associated with clinical outcome. To further characterize spatial and phenotypic changes at high resolution, artificial intelligence-assisted digital image analysis of hematoxylin and eosin sections and spatial transcriptomics (10x Genomics Visium) were performed on paired tissue samples. In baseline samples, cells residing in lymphoid aggregates and tertiary lymphoid structures (LAs/TLSs) were more abundant in stable disease than in relapse and exhibited gene expression programs associated with improved survival in UC. Most spatial features reflected shared pharmacodynamic changes between stable disease and relapse; however, carcinoma-endothelial adjacency was reduced significantly following treatment and differed between groups, accompanied by distinct transcriptional programs. Together, these findings indicate that atezolizumab induces localized immune and stromal remodeling within the tumor microenvironment, while non-response despite immune expansion is associated with persistent spatial immune exclusion and carcinoma-endothelial adjacency. Spatial and phenotypic biomarkers identified here may inform rational combination strategies for immune checkpoint inhibitor-refractory urothelial carcinoma.

PURPOSE: Newly-diagnosed glioblastoma (nGBM) is a devastating tumor with median survival of only 14-18 months despite aggressive standard of care (SOC). Dendritic cell (DC) homologous antigenic double-loading provides a powerful pattern-based signal that initiates cDC1-like skewing of monocytic precursors, inducing downstream development of CD8+ memory effectors. Here we report phase I results for DOC1021 (dubodencel), a novel DC vaccine regimen integrated with SOC. METHODS: In this dose-escalating study, DC prepared from mobilized peripheral blood were doubly loaded with autologous tumor lysate and amplified tumor mRNA and administered bilaterally near the deep cervical node chains in three biweekly courses given with weekly peg-IFN after conclusion of chemoradiation. Four dose levels from 3.5x106 to 3.6x107 total cells were tested. Patients with subtotal resection or tumor progression prior to vaccination were not excluded. RESULTS: Eighteen patients (median age 61 years (range 47-73), 94% MGMT unmethylated, 25% subtotal/partial resected) completed vaccination (16 nGBM, 2 recurrent) with no dose-limiting toxicities. Attributable AE were mostly mild and flu-like or injection-site reactions. Twelve-month OS among the newly-diagnosed cohort was 88% compared to an expected ~60% for SOC alone. Patients who received observation rather than reoperation in response to worsening MRI contrast-enhancement demonstrated gradual lesional resolution and improved OS. Immunophenotyping revealed post-vaccination elevations in CD4 and CD8 memory T-cells in peripheral blood, and spatial transcriptomic analysis revealed foci of activated inflammatory complexes at the primary tumor site. CONCLUSIONS: DOC1021 was safe, feasibly integrated within SOC, and associated with more favorable outcomes in this challenging patient population. Patients who received observation rather than reoperation for worsening MRI contrast-enhancement exhibited superior survival, suggesting an immune-reactive tumor microenvironment manifesting as pseudo-progression. These data supported initiation of a randomized Phase II trial (NCT06805305) for nGBM.

Background: Predicting whether a treatment will demonstrate meaningful clinical benefit before committing to a large-scale trial remains a major unmet need in oncology. Patient-derived organoids (PDOs) recapitulate individual tumor drug sensitivity, but have not been used to forecast population-level trial outcomes. We developed SCOPE (Screening-to-Clinical Outcome Prediction Engine), a platform that integrates PDO drug screening with clinical prognostic modeling to predict arm-level median progression-free survival (mPFS) and objective response rate (ORR) without access to any trial outcome data. Patients and methods: SCOPE was trained on 54 treatment lines from patients with metastatic colorectal cancer (mCRC, n=15) and metastatic pancreatic ductal adenocarcinoma (mPDAC, n=39) with matched clinical data and PDO drug screening across 9 compounds. A Clinical Score module captures baseline prognosis; a Drug Screen Score module quantifies treatment-specific organoid sensitivity. To predict trial outcomes, synthetic patient profiles are generated from published eligibility criteria and matched to a biobank of 81 PDO lines. Predictions were externally validated against 32 arms from 23 published trials, treatment ranking was assessed across 8 head-to-head comparisons, and prospective applicability was tested for daraxonrasib (RMC-6236), a novel pan-RAS inhibitor in mPDAC. Results: Predicted mPFS strongly agreed with published outcomes (R2=0.85, MAE=0.82 months; Pearson r=0.92, P<0.001), approaching the empirical concordance between two independently measured clinical endpoints (ORR vs. mPFS, R2=0.87). ORR prediction was similarly robust (R2=0.71, MAE=7.3 percentage points). Integrating organoid and clinical data significantly outperformed either alone (P=0.001). SCOPE correctly identified the superior arm in 7 of 8 head-to-head comparisons (88%, P<0.05). Applied to daraxonrasib prior to phase 3 data availability, the platform predicted superiority over standard chemotherapy in KRAS-mutant mPDAC, consistent with emerging clinical data. Conclusion: By combining functional organoid drug screening with clinical modeling, SCOPE generates calibrated efficacy predictions for both established regimens and novel agents without prior clinical data. This approach could support clinical trial design, treatment arm selection, and go/no-go decisions, offering a new tool to improve the efficiency of gastrointestinal cancer drug development.

Background: The tumor suppressor gene TP53 and the oncogene KRAS are among the most frequently altered core drivers in human malignancies. Although they cooperatively regulate critical biological processes, the prognostic impact of their co alterations remains poorly defined and exhibits striking inconsistency across different cancer types. Methods: We comprehensively analyzed genomic and clinical data from multi-cancer cohorts sourced from the cBioPortal database and The Cancer Genome Atlas (TCGA). Genetic alterations, including sequence variations and copy number alterations (CNAs), were classified for TP53 and KRAS. Patients were stratified into four subgroups based on individual or combined alteration status. Survival analyses were performed using Kaplan-Meier methods. Integrated multi-omics analyses were conducted to assess the relationship between genetic alterations and mRNA/protein expression, and to characterize co-occurring genetic events and their prognostic implications. Results: Patients harboring concurrent TP53 and KRAS alterations exhibited significantly shorter overall survival in pancreatic cancer, colorectal cancer, and ampullary carcinoma, but surprisingly demonstrated the longest survival in gastric cancer. Distinct KRAS mutation subtype distributions were observed across cancer types: G12D/G12V predominated in pancreatic and colorectal cancers, G12C in non small cell lung cancer, and G13D in gastric cancer, with copy number alterations representing a substantial proportion of KRAS alterations in gastric and lung cancers. Multi-omics analysis revealed a lack of concordance between genetic alterations and mRNA/protein expression, indicating that mutation status alone does not reliably reflect downstream molecular changes. Concurrent genetic events displayed striking cancer-type specificity: CDKN2A alterations frequently co-occurred with TP53/KRAS double alterations in pancreatic cancer and were associated with worse prognosis, whereas APC mutations co-occurred in colorectal cancer and correlated with improved survival. Integrated analysis further demonstrated that KRASaltered/TP53altered patients were highly enriched in pancreatic, colorectal, and lung cancers, each exhibiting unique background genomic landscapes. Conclusions: The prognostic significance of TP53 and KRAS alterations is profoundly cancer-type specific, driven by differences in mutation subtype distribution, copy number alteration patterns, co-occurring genetic events, and the discordance between genotype and functional expression. These findings challenge the simplistic view of dual-gene alterations as universal markers of poor prognosis and underscore the necessity of incorporating cancer-specific molecular contexts into prognostic models and precision oncology strategies.

Background. Pancreatic ductal adenocarcinoma (PDAC) has a five-year survival rate of approximately 12%, largely because it is typically diagnosed at an advanced stage. CT-based computational methods for early detection exist but rely on black-box deep learning or large texture feature sets without tissue-specific interpretability. Methods. We developed Virtual Spectral Decomposition (VSD), which applies six parameterized sigmoid functions S(HU) = 1/(1+exp(-alpha x (HU - mu))) to standard portal-venous CT, decomposing each pixel into tissue-specific response channels for fat (mu=-60), fluid (mu=10), parenchyma (mu=45), stroma (mu=75), vascular (mu=130), and calcification (mu=250). Dendritic Binary Gating identifies structural content per channel using morphological filtering, enabling co-firing analysis and lone firer identification. A 25-feature signature was extracted per patient. Three independent datasets were analyzed: NIH Pancreas-CT (n=78 healthy), Medical Segmentation Decathlon Task07 (n=281 PDAC, paired tumor/adjacent tissue), and CPTAC-PDA from The Cancer Imaging Archive (n=82, multi-institutional, with DICOM time point tags). The same six sigmoid parameters were used across all datasets without retraining. Results. VSD achieved AUC 0.943 for field effect detection (healthy vs cancer-adjacent parenchyma) and AUC 0.931 for patient-stratified tumor specification on MSD. On CPTAC-PDA, VSD achieved AUC 0.961 (6 features) and 0.979 (25 features) for distinguishing healthy from cancer-bearing pancreas on scans obtained prior to pathological diagnosis. All significant features replicated across datasets in the same direction: z_fat (d=-2.10, p=3.5e-27), z_fluid (d=-2.76, p=2.4e-38), fire_fat (d=+2.18, p=1.2e-28). Critically, VSD severity did not correlate with days-from-diagnosis (r=-0.008, p=0.944) across a range of day -1394 to day +249. Patient C3N-01375, scanned 3.8 years before pathological diagnosis, had VSD severity 1.87, well above the healthy mean of 0.94 +/- 0.33. The tissue transformation signature was temporally stable, indicating an early, persistent tissue state rather than a progressively worsening process. Conclusions. VSD with Dendritic Binary Gating detects a stable pancreatic tissue composition signature on standard CT that is present years before clinical diagnosis, validated across three independent datasets without parameter adjustment. The six sigmoid channels map to biologically meaningful tissue components through a fully transparent interpretability chain. The temporal stability of the signal implies a detection window of 3-7 years, consistent with known PanIN-3 microenvironment transformation timelines. VSD functions as a single-scan screening tool applicable to any abdominal CT performed during the pre-clinical window.

Background: Immune checkpoint inhibitors (ICIs) have revolutionized cancer therapy but can cause serious immune-related adverse events (irAEs), with pneumonitis (ICI-P) being among the most severe. Early identification of high-risk patients before ICI initiation is critical for closer monitoring, timely intervention, and improved outcomes. Purpose: To develop and validate a deep learning foundation model to predict ICI-P from baseline CT scans in patients with lung cancer. Methods: We designed the Checkpoint-Inhibitor Pneumonitis Hazard EstimatoR (CIPHER), a deep learning foundation model that combines contrastive learning with a transformer-based masked autoencoder to predict ICI-P from baseline CT scans in patients with lung cancer. Using self-supervised learning, CIPHER was pre-trained on 590,284 CT slices from 2,500 non-small cell lung cancer (NSCLC) patients to capture heterogeneous lung parenchymal patterns. After pre-training, the model was fine-tuned on an internal NSCLC cohort for ICI-P risk prediction, using images from 254 patients for model development and 93 patients for internal validation. We compared CIPHER with classical radiomic models and further evaluated it on an external NSCLC cohort of 116 patients. Results: In the internal immunotherapy cohort, CIPHER consistently distinguished patients at elevated risk of ICI-P from those without the event, with AUCs ranging from 0.77 to 0.85. In head-to-head benchmarking, CIPHER achieved an AUC of 0.83, outperforming the radiomic models. In the external validation cohort, CIPHER maintained strong performance (AUC = 0.83; balanced accuracy = 81.7%), exceeding the radiomic models (DeLong p = 0.0318) and demonstrating higher specificity without sacrificing sensitivity. By contrast, the radiomic model showed high sensitivity (85.0%) but markedly lower specificity (45.8%). Confusion matrix analysis confirmed the robust classification performance of CIPHER, correctly identifying 80 of 96 non-ICI-P cases and 16 of 20 ICI-P cases. Conclusions: We developed and externally validated CIPHER for predicting future risk of ICI-P from pre-treatment CT scans. With prospective validation, CIPHER may be incorporated into routine patient management to improve outcomes.

BackgroundTreatment optimization in HPV-associated oropharyngeal cancer (OPSCC) remains challenging, as recent de-escalation trials have shown limited success. Current patient selection strategies based on smoking history and TNM classification are insufficient, highlighting the need for robust, standardized prognostic biomarkers. We report the first validation of the Immunoscore (IS) for prognostic stratification in HPV-associated OPSCC. Patients and methodsWe analyzed 191 HPV-associated (p16+ and HPV DNA/RNA+) OPSCC patients from an international multicenter cohort (2015-2024), comprising a French monocentric retrospective training cohort (N = 48) and three validation cohorts: French monocentric retrospective (N = 48), French multicenter prospective (N = 50), and US multicenter retrospective (N = 45). IS is a standardized digital pathology assay quantifying CD3lJ and CD8lJ densities in tumor cores and invasive margins, with cut-offs defined in the training cohort and validated across cohorts. Associations with disease-free survival (DFS), time to recurrence (TTR) and overall survival (OS) were assessed, alongside 3RNA-seq and sequential immunofluorescence profiling of immune composition. ResultsMedian age 65; 80% male; 74% smokers; 66% T1-2; 82% N0-1 (AJCC8th). IS-High patients demonstrated superior 3-year DFS in the training and validation cohorts 1-3 (all log-rank P < 0.05). Multivariable analysis identified IS-Low as the strongest independent risk factor for DFS (HR 9.03; 95% CI: 4.02-20.31; P < 0.001). The model combining IS with clinical factors showed higher predictive accuracy for DFS (C-index 0.82) than clinical variables alone (0.7; P < 0.0001). Similar findings were observed for TTR and OS. IS-High tumors showed markedly higher enrichment of lymphoid and myeloid immune cell populations, contrasting with immune-poor signatures in IS-Low tumors. ConclusionsIS is a robust biomarker that outperforms standard clinical variables in both prognostic and predictive accuracy. The enriched cytotoxic immune infiltrate in IS-High tumors explains favorable outcomes and supports their suitability for treatment de-escalation. Prospective validation is warranted.

Background: Sezary syndrome (SS) is an aggressive leukemic variant of cutaneous T-cell lymphoma (CTCL) with distinct clinical and biological features compared to rarer entities such as primary cutaneous CD8+ aggressive epidermotropic cytotoxic T-cell lymphoma (PCAECTCL). Although recurrent genomic alterations in CTCL have been described, comparative analyses at the pathway level across biologically divergent subtypes remain limited. Here, we leveraged a conversational artificial intelligence (AI) platform for precision oncology to enable rapid, integrative, and hypothesis-driven interrogation of publicly available genomic datasets. Methods: We conducted a secondary analysis of somatic mutation and clinical data from the Columbia University CTCL cohort accessed via cBioPortal. Cases were stratified into SS (n=26) and PCAECTCL (n=13). High-confidence coding variants were curated and mapped to biologically relevant signaling pathways and functional gene categories implicated in CTCL pathogenesis. Pathway-level mutation frequencies were compared using Chi-square or Fisher's exact tests, with effect sizes quantified as odds ratios. Tumor mutational burden (TMB) was compared using the Wilcoxon rank-sum test. Subtype-specific co-mutation patterns were evaluated using pairwise association analyses and visualized through oncoplots and network heatmaps. Conversational AI agents, AI-HOPE, were used to iteratively refine cohort definitions, prioritize pathway-level signals, and contextualize findings. Results: TMB was comparable between SS and PCAECTCL (p = 0.96), indicating no significant difference in global mutational load. In contrast, pathway-centric analyses revealed marked qualitative differences. SS demonstrated enrichment of alterations in epigenetic regulators, tumor suppressor and cell-cycle control pathways, NFAT signaling, and DNA damage response mechanisms, consistent with transcriptional dysregulation and immune modulation. PCAECTCL exhibited relatively higher frequencies of alterations involving epigenetic regulators and MAPK pathway signaling, suggesting distinct oncogenic dependencies. Co-mutation analysis revealed a more constrained and focused interaction landscape in SS, whereas PCAECTCL displayed broader and more heterogeneous co-mutation networks, indicative of divergent evolutionary trajectories. Notably, ERBB2 mutations were significantly enriched between subtypes (p = 0.031), highlighting a potential subtype-specific therapeutic vulnerability. Conclusions: This study demonstrates that SS is distinguished from PCAECTCL not by increased mutational burden but by distinct pathway-level architectures, particularly involving epigenetic regulation, immune signaling, and transcriptional control. These findings generate biologically grounded, testable hypotheses for subtype-specific therapeutic targeting and underscore the value of conversational AI as a scalable framework for accelerating discovery in translational cancer genomics.

Background: Current deep learning models in computational pathology, radiology, and digital pathology produce opaque predictions that lack the explainable artificial intelligence (xAI) capabilities required for clinical adoption. Despite achieving radiologist-level performance in tasks from whole-slide image (WSI) classification to mammographic screening, these models function as black boxes: clinicians cannot trace predictions to specific biological features, verify outputs against established morphological criteria, or integrate AI reasoning into precision oncology workflows and tumor board decision-making. Methods: We present Virtual Spectral Decomposition (VSD), a modality-agnostic, interpretable-by-design framework that decomposes medical images into six biologically interpretable tissue composition channels using sigmoid threshold functions - the same mathematical structure as CT windowing. Unlike post-hoc xAI methods (Grad-CAM, SHAP, LIME) applied to black-box deep learning models, VSD channels have pre-defined biological meanings derived from tissue physics, providing inherent explainability without sacrificing quantitative rigor. For whole-slide image (WSI) analysis in digital pathology, we introduce the dendritic tile selection algorithm, a biologically-inspired hierarchical architecture achieving 70-80% computational reduction while preferentially sampling the tumor immune microenvironment. VSD is validated across three cancer types and imaging modalities: pancreatic ductal adenocarcinoma (PDAC) on CT imaging, lung adenocarcinoma (LUAD) on H&E-stained pathology slides using TCGA data, and breast cancer on screening mammography. Composition entropy of the six-channel vector is computed as a visual Biological Entropy Index (vBEI) - an imaging biomarker quantifying the diversity of active biological defense systems. Results: In pancreatic cancer, the fat-to-stroma ratio (a novel CT-derived radiomics biomarker) declines from >5.0 (normal) to <0.5 (advanced PDAC), enabling early detection of desmoplastic invasion before mass formation on standard imaging. In lung cancer, composition entropy from H&E whole-slide images correlates with tumor immune microenvironment markers from RNA-seq (CD3: rho=+0.57, p=0.009; CD8: rho=+0.54, p=0.015; PD-1: rho=+0.54, p=0.013) and predicts overall survival (low entropy immune-desert phenotype: 71% mortality vs 29%, p=0.032; n=20 TCGA-LUAD), providing immune phenotyping for checkpoint immunotherapy patient selection from a $5 H&E slide without molecular assays. In breast cancer, each lesion type produces a characteristic six-channel fingerprint functioning as an interpretable computer-aided diagnosis (CAD) system for quantitative BI-RADS assessment and subtype classification (IDC vs ILC vs DCIS vs IBC). A five-level xAI audit trail provides complete traceability from clinical decision support output to specific biological structures visible on the original images. Conclusion: VSD establishes a unified, interpretable-by-design mathematical framework for explainable tissue composition analysis across imaging modalities and cancer types. Unlike black-box deep learning and post-hoc xAI approaches, VSD provides inherently interpretable, clinically verifiable cancer detection and immune phenotyping from standard clinical imaging at existing costs - without requiring foundation model infrastructure, specialized hardware, or molecular assays. The open-source pipeline (Google Colab, Supplementary Material) enables immediate reproducibility and extension to additional cancer types across the pan-cancer TCGA atlas.

Soft tissue leiomyosarcoma (STLMS) is an aggressive malignancy for which robust molecular subclassification and mechanism-based therapeutic strategies still remain limited. We performed integrative proteogenomic analyses of primary and metastatic STLMS to define subtype-associated molecular programs. Joint analysis of the proteome and phosphoproteome identified 3 biologically distinct subtypes. P1 was characterized by relative genomic stability, low proliferative activity, and enrichment of FGFR2- and PDK-associated signaling. In contrast, P2 and P3 showed greater chromosomal instability and more aggressive clinical behavior, but with distinct molecular features. Notably, P2 was associated with inflammatory and RTK-RAS pathway programs, activation of CDK-AURKA/B-mTOR-ERK kinase networks, IGF1R/PDGFRA alterations, and the poorest outcomes. On the other hand, P3 showed strong cell cycle and DNA repair programs, elevated NCOR1 expression, and increased expression of nonhomologous end joining components, including PARP1. Homologous recombination deficiency analyses distinguished HRD-low P1 from HRD-high P2/P3, and paired analyses suggested increased HRD-related features in metastatic lesions within P3. Immune profiling identified an immune-hot yet potentially suppressive state in P2, marked by higher LGALS9 expression and M2-like macrophage infiltration. To support clinical translation, we developed a tissue microarray-based immunohistochemical classifier that enabled surrogate assignment of proteome-defined subtypes in an independent cohort and showed recurrence-free survival differences across inferred subtypes. These findings together establish a proteogenomic framework for STLMS heterogeneity and nominate subtype-associated biological vulnerabilities for future translational and clinical investigation.

BackgroundNeoadjuvant chemoimmunotherapy (nCIT) has become a standard treatment for locally advanced resectable non-small cell lung cancer (NSCLC), yet the spatial biology underlying treatment resistance remains poorly understood. We used spatial transcriptomics to define the microenvironmental architecture of residual cancers in patients who did not achieve major pathologic response (non-MPR) compared with those who did (MPR). MethodsSpatial transcriptomics was performed on 10 formalin-fixed paraffin-embedded (FFPE) tumor blocks (5 MPR, 5 non-MPR) obtained from 8 patients treated with nCIT. A deep learning algorithm was applied to detect viable residual cancer spots from treatment-induced fibrosis and necrosis. Spatial deconvolution, distance modeling, ligand-receptor analysis, and functional pathway scoring were integrated to characterize niche-specific programs. ResultsMPR cancer core displayed an immune-permissive remodeling environment with deep infiltration of cytotoxic CD8+ T cells, mature dendritic cells (LAMP3+, CCR7+), and active efferocytosis signaling (APOE-TREM2), alongside robust MHC class II expression. Non-MPR cancer core, by contrast, exhibited spatial immune exclusion: a dense fibroblast barrier reinforced by TIMP1-CD63 signaling and Treg-enriched boundaries physically restricted effector T cell access to the cancer core. Residual cancer cells in non-MPR samples maintained active cell cycling and independently upregulated cytochrome P450-mediated drug detoxification and DNA damage response pathways without inducing MHC class II expression -- effectively decoupling intrinsic survival from immune recognition. The non-MPR core also showed a hyper-metabolic profile, including elevated glutathione metabolism consistent with antioxidant buffering against chemotherapy-induced oxidative stress. TROP2 was broadly expressed across the non-MPR cancer core and co-localized with DNA damage response and nuclear factor erythroid 2-related factor 2 resistance signatures. ConclusionsResidual cancer cores in non-MPR tumors appear to represent evolved resistant niches sustained by structural immune exclusion, metabolic rewiring, and DNA repair proficiency. These findings highlight the spatial co-localization of epithelial anchors, such as TROP2, with intrinsic resistance pathways, providing a structural rationale for developing novel precision therapeutic strategies to bypass stromal barriers and overcome the cancer cores intrinsic repair capacity.

Keratin-positive giant cell-rich tumor (KPGCT) is a newly described bone and soft tissue tumor. The tumor is characterized by scattered keratin-positive cells and the presence of HMGA2::NCOR2 fusions. It is not known if the HMGA2::NCOR2 fusion is located in the keratin-positive cells, and little is known about how KPGCT develops. KPGCT shares some histologic features with tenosynovial giant cell tumor (TGCT), a soft tissue tumor with CSF1 rearrangements. Single-nuclei RNA sequencing (snRNA-seq) and Xenium spatial transcriptomics were used to elucidate the mechanisms driving KPGCT and compare KPGCT to TGCT. We show that the neoplastic cells in KPGCT constitute only a minority of cells in the tumor, and that they co-express keratin, HMGA2 and CSF1. The neoplastic cells in KPGCT express no synovial markers, confirming KPGCT as a distinct entity, separate from TGCT. The bulk of the tumor consists of CSF1R-expressing macrophages and osteoclast-like giant cells, suggesting an important role for CSF1-CSF1R signaling. In addition, we find that the cells with the HMGA2 translocation show activation of the hippo signaling pathway, which is known to regulate CSF1 expression. We show that the CSF1-CSF1R axis, possibly regulated through the hippo signaling pathway, plays an important role in KPGCT. This axis likely stimulates the migration and proliferation of macrophages, which form the majority of cells in the tumor, as well as their differentiation into osteoclasts-like giant cells. These results provide a rationale for the use of CSF1R inhibitors, which have already shown efficacy in TGCT, as a therapy for KPGCT. SignificanceKeratin-positive giant cell-rich tumor (KPGCT) is a rare, newly described soft tissue and bone tumor. By examining this tumour on a single-cell level, we confirm the identity of the neoplastic cells on a molecular level, showing these form a minority of cells in the tumor. We show that activation of the hippo pathway in the neoplastic cells is a likely driver of tumorigenesis. Additionally, we show the neoplastic cells produce large amounts of CSF1, attracting the macrophages that form the majority of cells in the tumor. This finding gives supporting evidence for anecdotal reports of response to CSF1 inhibitor therapy. Finally, we identify key differences between KPGCT and tenosynovial giant cell tumor, a tumor that shares histological features with KPGCT.

Purpose: To test whether histology-derived gene-expression signatures from routine hematoxylin and eosin slides are prognostic for recurrence and predictive of chemotherapy benefit in early breast cancer. Methods: We conducted a multi-cohort study including CALGB 9344 (anthracycline +/- paclitaxel), CALGB 9741 (standard vs dose-dense chemotherapy), a pooled Chicago real-world cohort, and the American Cancer Society (ACS) Cancer Prevention Studies-II and -3. Whole-slide images were processed with a previously described pipeline to generate 61 histology-derived signatures per patient. The primary endpoint was distant recurrence-free interval (DRFI), except in ACS, where breast cancer-specific survival was used. Secondary endpoints include distant recurrence-free survival (DRFS) and overall survival. The most prognostic signature in CALGB 9344, selected by Harrell's C-index, was evaluated in additional cohorts. Signature-treatment interaction was assessed by likelihood-ratio tests. Multivariable Cox models incorporating age, tumor size, nodal status, estrogen/progesterone receptor status, and signature were fit in CALGB 9344 to improve risk stratification. Results: A total of 7,170 patients were included across four cohorts. The top histology-derived signature in CALGB 9344 showed strong prognostic performance for 5-year DRFI (C-index 0.63) and performed well across validation cohorts (C-index 0.60, 0.70, and 0.62 in CALGB 9741, Chicago, and ACS, respectively). The strongest predictive signal for treatment benefit was observed for DRFS. High-risk cases identified by the signature demonstrated greater benefit from taxane in CALGB 9344 (adjusted hazard ratio [aHR] 0.76 for DRFS, 95% CI 0.66-0.88; interaction p=0.028), from dose-dense chemotherapy in CALGB 9741 (aHR 0.69, 95% CI 0.56-0.85; interaction p=0.039), and differential chemotherapy benefit in the Chicago cohort (aHR 0.84, 95% CI 0.59-1.21; interaction p=0.009). Combined clinical-histology models improved risk stratification and identified low-risk groups with a 2%-10% risk of distant recurrence or breast cancer death. Conclusion: Histology-derived signatures from H&E images are broadly prognostic and, unlike clinical factors, may predict chemotherapy benefit.

KIF18A inhibition has emerged as a therapeutic strategy for chromosomally unstable cancers, but clinical development is limited by the absence of a deployable predictive biomarker. Here we identify strong, diffuse p16INK4a expression, a well-established surrogate marker of Rb-pathway inactivation, as a predictive biomarker of response to KIF18A inhibition, and show that Rb-pathway inactivation marks a biologically distinct subset of cancers sensitive to this therapeutic approach. In sensitive models, low Rb activity is associated with robust spindle assembly checkpoint signaling and prolonged mitotic arrest following KIF18A inhibition. Weakening the spindle assembly checkpoint in this context is sufficient to confer resistance. Across three independent pan-cancer sensitivity datasets generated with distinct KIF18A inhibitors, Rb-pathway altered models were significantly more sensitive than histology-matched Rb-intact comparators, with the strongest association observed in cancers harboring direct RB1 loss or inactivating mutation. Guided by this mechanism, we retrospectively analyzed p16INK4a expression by immunohistochemistry (IHC) in pre-treatment tumor biopsies from 79 heavily pre-treated high-grade serous ovarian cancer patients across three dose-escalation or expansion cohorts and treated with two different KIF18A inhibitors (sovilnesib and VLS-1488) sharing a common mechanism of action. p16INK4a-high tumors showed substantially higher objective response rates than their p16INK4a-low counterparts (36.0% versus 2.2%; P = 0.0002) and markedly longer progression-free survival (median 24.3 versus 7.9 weeks; hazard ratio, 0.16; P < 0.0001). These findings establish p16INK4a as a mechanistically-based, clinically implementable biomarker of clinical response to KIF18A inhibition that is poised to support pan-cancer development of KIF18A inhibitors guided by Rb-pathway inactivation.

We studied the effect of stimulating innate immune function in tumor-associated myeloid cells (TAMs) in medulloblastoma (MB) and diffuse midline glioma (DMG), using a polyoxazoline nanoparticle formulation of the TLR7/8 agonist resiquimod (ResiPOx). Children with MB and DMG need novel therapeutic strategies to improve outcomes and reduce recurrence. We investigated the effect of systemically administered ResiPOx on TAMs in MB and DMG using endogenous MB and DMG models in immune-competent mice and identified multiple mechanisms of anti-tumor effect. We packaged resiquimod into polyoxazoline micelles to generate ResiPOx. We studied ResiPOx efficacy as a single agent or paired with radiation therapy (RT). We determined ResiPOx pharmacokinetics (PK) using tritium-labeled resiquimod and mass spectroscopy imaging (MSI). We determined ResiPOx pharmacodynamics (PD) using flow cytometry immunohistochemistry, bulk and single-cell RNA-seq and immunoblotting. We then studied ResiPOx safety and PD in a non-human primate model using rhesus macaques. ResiPOx formulation improved the blood-brain barrier penetration and anti-tumor efficacy of resiquimod. ResiPOx treatment extended progression-free survival (PFS) in mice with MB and DMG. In both tumor types, ResiPOx expanded TAM populations and reprogrammed TAMs toward anti-tumoral states, characterized by activation of IFN{beta} and extrinsic apoptosis pathway signaling, antigen presentation, and T cell activation signatures. In rhesus macaques, systemic ResiPOx administration was well tolerated and induced brain transcriptional responses that resembled ResiPOx responses in DMG and MB mouse models, indicating common effects across species from mice to non-human primates, and highlighting potential for similar effects in patients. ResiPOx is a brain-penetrant immunomodulatory therapeutic that reshapes the immune-privileged brain tumor microenvironment. Systemic administration activates myeloid-driven anti-tumoral immunity mediated by microglial and macrophage TAMs, and improves survival in preclinical models of DMG and MB.

Mechanistic modeling has long been used as a tool to describe the dynamics of biological systems, especially cancer in response to treatment. Their key advantage lies in interpretability of relationships between input parameters and outcomes of interest. In contrast, machine learning techniques offer strong prediction performance, especially for high dimensional datasets that are common in oncology. Here, we employ a Mechanstic Learning framework that combines the advantages of both approaches by training machine learning models on mechanistic parameters inferred from clinical patient data. The mechanistic model (a Markov chain model) contains sixteen parameters that describe the rate of cell fate transitions that occur in patients with B-cell precursor acute lymphoblastic leukemia. The machine learning (a ridge logistic regression model) is trained on these parameters to predict two clinically-relevant features: BCR::ABL1 fusion gene status (positive or negative) and minimal residual disease status (positive or negative) post-induction chemotherapy. Model training is done in an iterative fashion to assess which (and how many) parameters are critical to maintain high predictive performance. Using machine learning models trained on the clinical flow-cytometry data, we find that the stem-like cell state alone is the most predictive feature for both BCR::ABL1-positive and MRD-positive disease, with combination scores (defined as the average of accuracy, balanced accuracy, and area under the curve) of 0.80 and 0.67, respectively. By comparison, mechanistic learning achieves comparable or improved combination scores for BCR::ABL1-positive and MRD-positive disease, with scores of 0.81 and 0.71, respectively, using only de-differentiation for BCR::ABL1 and primitive-state persistence together with differentiation-directed exit for MRD. Thus, the mechanistic-learning approach not only preserves predictive performance, but also provides a biological hypothesis for why stemness is predictive of these clinically relevant outcomes.

Background: Stomach adenocarcinoma is driven by heterogeneity, limiting therapeutic success. Although ROS acts as a continuous redox rheostat for tumor evolution, it is categorized based on binary models that are masked by tumor-microenvironment (TME) confounders. Here, we have defined a continuous, TME-independent ROS axis to help identify intrinsic vulnerabilities and improve patient stratification. Methods: Non-negative matrix factorization (NMF) defined a ROS-Axis in TCGA-STAD which was validated in ACRG Cohort. Multivariate regression model isolated intrinsic signatures via residual ROS scores by adjusting for TME confounders. Survival was assessed using Cox hazard models. Drug sensitivities were mapped using GDSC2/ElasticNet modeling with cross-cohort replication. Results: Our results define a reproducible ROS gradient, driven by effectors like NQO1 and SOD1, characterizing ROS-high tumors as proliferative, epithelial and immune -cold. High residual ROS score was associated with an improved prognosis, regardless of TNM stage and age. Pharmacogenomic mapping revealed an overlapping sensitivity to mTOR inhibitors in ROS-high gastric cancer tumors which persisted after TME confounder adjustment. Conclusion: The continuous ROS axis provides a functional readout of metabolic dependency that refines traditional anatomical staging. By identifying mTOR dependent cold tumors, our framework offers a precision strategy for immunotherapy-resistant patients like those affected by microsatellite-stable gastric cancer.

Despite considerable advances, the emergence of treatment resistance to tyrosine kinase inhibitors (TKIs) therapy remains a significant challenge in chronic myeloid leukemia (CML). Here, we report the first clinical case of resistance to combined ponatinib and asciminib therapy in a CML patient who relapsed with B lymphoblastic blast crisis. While at presentation the patient harbored the canonical e13a2 BCR::ABL1 fusion, at relapse his disease harbored the T315I mutation together with a novel e6a3 BCR::ABL1 fusion, arisen by internal deletion in the original translocated allele. Structural modeling and biochemical analyses demonstrated that deletion of exon 2-encoded residues of ABL1 destabilizes the autoinhibited conformation, resulting in a hyperactive kinase with increased propensity for B-cell differentiation. Functional studies revealed that both BCR::ABL1e6a3 and BCR::ABL1e6a3/T315I conferred resistance to ponatinib and asciminib, alone or in combination. BCR::ABL1e6a3 demonstrated enhanced sensitivity to active-state selective inhibitors dasatinib and bosutinib, whereas BCR::ABL1e6a3/T315I remained resistant. Combined drug sensitivity assays showed that axitinib restored inhibitory activity when combined with ponatinib or asciminib. Strikingly, a combination of axitinib and asciminib with low dose ponatinib fully suppressed enzymatic activity of BCR::ABL1e6a3/T315I and cellular proliferation. These data show that treatment with asciminib and ponatinib can select for mutations with notably elevated enzymatic activity, effectively targeted by an axitinib-based triple combination. These data highlight the remarkable mutability of the BCR::ABL1 kinase, including through novel isoforms and provides a strong rationale for the clinical assessment of a triple inhibitor combination as a strategy to overcome resistance to dual ponatinib and asciminib therapy.

Protein arginine methyltransferase 5 (PRMT5) is a synthetic lethal target in methylthioadenosine phosphorylase-deleted (MTAP-null) cancers. Second-generation MTA-cooperative PRMT5 inhibitors preferentially target MTAP-null cells while largely sparing MTAP-wildtype (MTAP-WT) cells, thereby improving tumor selectivity over first-generation PRMT5 inhibitors. Despite encouraging efficacy and safety signals in early clinical studies, the modest objective response rates (ORRs) observed with these inhibitors suggest that intrinsic or acquired resistance may limit their clinical benefit. Here, we investigated mechanisms of acquired resistance to the MTA-cooperative PRMT5 inhibitor BMS-986504/MRTX1719 in MTAP-null non-small cell lung cancer (NSCLC) cells and sought to identify therapeutic vulnerabilities that emerge upon resistance. Using multiple in vitro-derived resistant models, we found that acquired resistance was not fully explained by alterations in PRMT5 activity or reduced MTA levels. Instead, resistance was associated with collateral sensitivity to MEK inhibition and enrichment of MAPK-related transcriptional programs. Together, these findings identify MEK inhibition as an actionable collateral vulnerability in MTAP-null NSCLC cells that acquire resistance to PRMT5 inhibition.